A Proposed High-Freqnency, 
Negative-Resistance Diode 

By W. T. READ, JR. 

(Miiriusrvipt Iteceived Miii-ch 2, 1957) 

This paper describes and analyzes a proposed semiconductor diode de- 
signed to operate as an oscillator when mounted in a suitable microwave 
cavity. The frequency would be in the range extending from 1 to 50 kmc. The 
negative Q may be as low as 10 and the efficiency as high as 30 per cent. 

The diode is biased in reverse so as to establish a depletion, or space- 
charge, layer affixed width in a relatively high resistance region, boundedby 
very low resistance end regions. The electric field has a maximum- at one 
edge of the space-charge region, where hole-electron pairs are generated by 
internal secondary emission, or avalanche. The holes {or electrons) travel 
across the space-charge layer nnth constant velocity, ih\is producing a cur- 
rent through the diode. Because of the build-up time of the avalanche, and 
the transit time of the holes across the depletion layer, the alternating current 
is delayed by approximately one-half cycle relative to the ac voltage. TMis, 
power is delivered to the ac signal. When the diode is mounted in an induc- 
tive microwave cavity tuned to the capacity of the diode, an oscillation will 
build up. It appears possible to obtain over 20 watts of ac power in con- 
timious operation at 6 kmc. 

I. DESCniPTION 

This paper discusses a proposed oscillator consisting of a semicon- 
ductor diode biased in reverse and mounted in a microwave cavity. The 
impedances of the cavity is mainly inductive and is matched to the 
mainly capa(atati\-e impedance of the diode so as to form a resonant 
system. We shall show that the diode can have a negati^'e ac resistance 
so that it delivers power from the dc bias to the oscillation. The negative 
Q may be as low as 10 and the efficiency as high as 30 per cent. 

The principle of operation is as follows: a reverse bias is applied to 
establish a space-charge, or depletion, layer of fixed width in a relatively 
weakly doped region bounded by highly doped end regions. A possible 
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stmctuic is the n"''-p-i-p''' structure, where the + denotes high doping 
and i means intrinsic. This structure is shown in Fig. l(:i). The field 
distribution under reverse bias is shown in Fig. 1(b). The voltage is al- 
ways well above the punch-through voltage, so that the space-<^hiirge 
region always extends from the n''"-p junction through the p and i (in- 
trinsic) regions to the i-p"^ junction. The fixed charges in the various re- 
gions are shown in Fig. 1(b). A positive charge gives a rising field in 
going from left to right. A positive field makes holes move to the right. 
The maximum field, which occurs at the n''"-p junction, is of the order of 
several hundred iviiovolts per cm, so that hole-electron pairs are gener- 
ated by internal secondary emission (also called multiplication or ava- 
lanche). The electrons go immediately into the n"^ region. The holes 
move to the right across the space-charge region. The field throughout 
the space-charge region is in the range (above about 5 kilovolts per cm) 
where carriers move with constant velocity independent of the field. ]''or 
practical purposes we can forget about the electrons in the space-eliargo 
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Fig. 1 — The structure (a) and field distribution (b) under reverse bias. 
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region and consider that current is generated at the n -p junction, so 
that only the holes move in the space-charge region. Thus, the physical 
picture is as follows: a current of holes h{t) is generated at the n'''-p 
junction. The holes move to the right and traverse the space-charge 
region, moving with constant velocity, v. In siUcon v is about 10^ cm sec" . 
The ti'ansit time of a hole across the space-charge region is r = W/v, 
where W is the width of the space-charge region. Throughout the dis- 
cussion we shall illustrate various quantities by giving numerical values 
for the case W = 10'^ cm. Thus t would be 10~' sec. 

The holes moving across the space-charge region produce a current 
hit) in the external circuit. It can be shown that h is equal to the aver- 
age current in the space-charge region. Since velocity is constant, /« is 
simply v/W = 1/t times the total charge of the movuig holes. (We shall 
assume throughout that all quantities refer to iniit area of junction.) 
Suppose, for example, that a pulse of holes of charge 5Q is suddenly gen- 
erated at the n'''-p junction. Immediately, a constant current Ic ^ 5Q/t 
begins to flow in the external circuit, and continues to flow during the 
time, T, that the holes are moving across the space-charge region from 
the n'^-p junction to the p"^ region. Thus, on the average, the external 
current Ie{t) due to the moving holes is delayed by t/2 relative to the 
current Io{t) generated at the n'^-p junction. We shall show (in the dis- 
cussion of multiplication) that the current Io{t) is delayed by one-quarter 
of a cycle relative to the ac voltage. Thus, to get a total delay of one-half 
cycle, we want the delay t/2 due to transit time to be one-quarter of a 
cycle. The cavity shoidd therefore be tuned to give a resonant frequency 
(2t)-\ 

The conductive current Te{t), which arises from carriers moving 
through the space-charge layer, should l)e distinguished from the dis- 
placement, or capacitative, current h , which charges and discharges 
the diode regarded as a capacitor. This current, h , supplies the vari- 
ation in charge at the edges of the space-charge region, where the field 
changes abruptly and, for practical purposes, can be considered discon- 
tinuous. Since Ic is 00° out of phase with the voltage it contributes 
nothing to the power. 

i.L Multiplication 

Carriers moving in the high field near the n"'"-p junction acquire 
enough energy to knock valence electrons into the conduction band, 
thus producing hole-electron pairs. The rate of pair production, or mul- 
tiplication, is a sensitive nonlinear fun(^tion of the field. Byproper doping, 
the field can be given a relatively sharp peak so that multiplication is 
confined to a very narrow region at the n"''-p junction. The multipUcation 
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rate can be regarded as a function of the peak field, Eo . Ato. critical field, 
Ec , breakdown will occur ; that is, any current will be self-sustaining ; 
every pair produced will on the average produce one other pair. We shall 
neglect the thermally generated i-everse saturation current, which will 
be small compared to the generated current Io{t). Then, when the field 
is above Ea , the current I<i(t) will be more than self-sustaining and will 
build up. When the field is below the critical breakdown field Ec , the 
current is less than self-sustaining and wiU die down. 

In operation, the diode is biased so that the peak field is above E^ 
during the positive half of the voltage cycle and below Ec during the 
negative half. Hence the current I(t) builds up during the positive half 
and dies down during the negative half. Therefore Io{i) reaches its maxi- 
mum in the middle of the voltage cj'-cle, or one-quarter of a cycle later 
than the voltage. 

We have assumed so far that the field varies in phase with the voltage. 
This will be a good approximation if the carrier space-charge can be 
neglected. This can be seen from Fig. 1(b). If the voltage is above the 
punch-through voltage, an increase in voltage simply raises the whole 
field distribution throughout the depletion layer. That is, additional 
charges simply appear at the edges of the space-charge layer. Hence the 
field at each point varies in phase with the voltage. The space charge of 
the carriers will, however, affect the shape of field distribution m the 
space-charge region. When the current becomes too large, the carrier 
space-charge cannot be neglected. As we shall see, the effect of carrier 
space-charge is to reduce the delay between the voltage and the current 
generated, 7o . In practice this limits the dc bias current. 

When the dc bias is small enough so that carrier space charge can be 
neglected, the operation can be summarized as follows: the peak field 
varies in phase with the voltage and generates (at one edge of the de- 
pletion layer) a current delayed by 90°. This current gives rise to a cur- 
rent through the diode delayed by 90° relative to the current generated 
and, therefore, by 180° relative to the voltage. 

1.2 Diode and Cavity 

The capacity of the diode and the inductance of the cavity determine 
the frequency, /, of oscillation. We have seen that the optimum fre- 
quency, /, is such that the transit time, t, is one-half cycle. Hence the 
angidar frequency, oj, should be 

For W = 10^^ cm, the frequency, /, would be 5,000 me. 
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We have seen that the current hito the diode is /« + /c , where I^ is ei 
(conductive current due to holes moving in the space-charge region and 
7c is a displacement current required to charge the diode regarded as a 
capacitor. Thus, the diode acts like a capacity and negative resistance in 
parallel. 

The (capacity is related to the width, W, of the depletion layer by 



a = 



1.06 X 10" 

W 



(2) 



for Si, 



where C is in farads per cm^ and W in cm. The negative conductance of 
the diode will be small compared to ojC, so the admittance is mainly that 
of the capacity. Thus for W = 10^^ cm the admittance will be about 
30 mhos per cm'. 

To make a resonant system we mount the diode in a microwave cavity 
having a mainly inductive impedance. The inductance, L, of the cavity 
is chosen so that the cavity and diode together have a resonant frequency 

Fig. 2 shows a possible ca\dty with the diode mounted in it. The diode 
should have narrow end regions so that the depletion layer is as close as 
possible to the metal base and center post of the cavity. Then the heat 
generated will flow away rapidly. The cavity is tuned by adjusting the 
vertical position of the plunger forming the top of the cavity. There 
must he some break in the cavity so that a dc bias can be applied. This 
is accomphshed in Fig. 2 by insulating the center lead. The break could 
also be made by insulating the walls from the base. On the right is an 
outlet through which the ac power is tapped. 
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Fig. 2 — Diode mounted in microwave cavity. 
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The impedance, wL, of the cavity must be equal to the impedance of 
the diode. If the diode is circular with radius R, then the junction area 
is ttR'^ and the impedance, {tR^uC)~^. Thus the area of the diode is 
limited by how small the impedance of the cavity can be made. Setting 
oiL = {wR'^wCy^ and using (1) and (2) we have 

'RY 4 



W / irKV(j}L ' , ^ 



0.5 X 10' 



for Si, 



where wL is in ohms. It would be difficult to make wL less than about 
10 ohms. Hence R could not be more than about 301^. 

1.3 Stable Operation 

We have shown that the diode will deliver power, so an oscillation will 
build up. In Section III we shall calculate exactly the Q for the small sig- 
nal case w^here everything is linear. For reasonable values of the direct cur- 
rent bias, Id , the Q can be made negative and as low in magnitude as 
desired. As the oscillation builds up, the behavior becomes highly non- 
hnear. In the practical range of operation the field will vary by a frac- 
tion — probably less than 20 per cent — of its average value but the 
current will vary by orders of magnitude. The current le approaches a 
square wave, being negligibly small during the positive half of the ac 
voltage cycle and almost constant during the negative half cycle. Since 
the direct current Id is the average conductive current, it follows that 
the amplitude of variation of le is approximately equal to Id . If Va is 
the ac voltage amplitude, the ac power delivered Pr is found to be 

P. = - hV. (4) 

TT 

per unit area of junction. Thus, if the dc bias is apphed by a constant 
direct current generator, the power delivered is proportional to Va ■ 
The energy stored in the capacity is proportional to Vj. Hence, — Q is 
proportional to Va ■ Iii other w^ords, if the amplitude increases, the 
stored energy, or energy of oscillation, increases faster than the energy 
delivered per cycle. This is the condition required for a stable oscillation 
to be possible. At the stable operating pomt, — Q of the diode is equal 
to the effective Q of the cavity. If the amphtudc increases, the energy 
dehvered by the diode mcreases less than the energy lost to the cavity. 
Thus, the amphtude decreases. In the same way a decrease in amplitude 
is also self-correcting. 

The total power output is Pr times the area irR^ and is therefore pro- 
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portional to R^IdV„ ■ The maximum R is given by (3). The maximum 
allowable la is limited by the fact that carrier space-charge cannot be 
too large, as we now sho\\'. 

1.4 Space Charge of the Carriers 

The holes moving across the space-charge region will affect the field 
distribution. As shown in Fig. 1 (b) a positive charge gives a rise in field 
(in going from left to right) and a negati^-c charge, a decrease in field. 
Thus, the effect of the holes will be to oppose the fixed negative charge 
in the depletion layer, and hence to flatten the field distribution shown 
in h'ig. 1(b)- Therefore, for a given voltage, the effect of the holes is to 
reduce the peak field at the n^-p junction. The multiphcation rate in- 
creases with the peak field. Thus, the current generated has a space 
charge that tends to reduce the rate of current generation. That is, the 
current tends to shut itself off. Instead of building up throughout the 
positive half of the voltage cycle, h(t) builds up until the carrier space- 
charge has reduced the peak field below the sustaining field E^ . Then the 
current decreases. Thus, /o reaches its peak hejore the middle of the 
cycle. This reduces the delay and hence the power. Increasing the cur- 
rent, therefore, increases the ac power only up to a certain point where 
carrier space-charge begins to spoil the phase relations. 

We have seen that in the practical range of operation the current /. 
varies by approximately its dc value Id ■ A current h corresponds to a 
total carrier charge rid in the space-charge layer (since I, is the average 
current in the space-charge layer and all carriers travel with velocity 
y = W/t). If the carrier space-charge is to be neglected it must be small 
compared to the charge CW that produces the voltage variation. Thus 
we want rid « CV^ ■ Actually, as we shall show in Section III, it is suffi- 
cient to take hr = CVJ'2. Since the period of oscillation is 2r, this be- 
comes 

Id = ~CV.. (5) 

If Id is no greater than this the current and voltage will be I'oughly 
ISO" out of phase and we can use (4) for the power per unit area. For 
larger Id the increase in power would be more than offset by the effect 
of carrier space-charge on the phase relations. 

1.5 Power Output 

Combining (3) with (o) we ha\-c, for the total ac power output 
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As we have seen, the impedance wL of the cavity could be as low as 10 
ohms. Now consider the maximum Va ■ This is limited by the dc voltage 
Vd . It is seen from Fig. 1(b) that if Va is too large the field in the in- 
trinsic region will be reduced to zero during the negative voltage cycle. 
By careful doping, the voltage at punch-through can be made small 
compared to the dc bias at breakdown. Even allowing for the reduction 
in field due to the carrier space-charge, we will be safe with Va = Vd/2. 
This gives an efficiency of about 30 per cent. From (6) the power be- 
comes 

Power Output - ttK'P. - (^] i ■ (7) 

\2t/ oiL 

The dc voltage is hmited by the following considerations: (a) In the 
negative half cycle the field must not fall below about 5 kilovolts per 
cm. Otherwise carrier velocity will depend on field; so the carriers will 
slow down during the negative half of the ac voltage cycle and thus re- 
duce the power. This can be avoided if Vd is at least lO^TF volts, where W 
is in cm. (b) To localize the multiphcation the field must remain well 
below Ec except near the n''"-p junction. This will be so if Fd = 2Va is 
less than about {2/^)ErW. Then the maximum field in the intrinsic 
region will always be less than 0.6 Ec . Thus we have 10*1F < Vd < 
0.4 EM. 

For TF as large as 10~^ cm the critical field will be about 350 kilovolts 
per cm. I'or IF as low as 10"'* cm the multiplication would have to be 
confined to a region no more than 10~ cm in width. This would require 
a somewhat steeper field gradient and a higher critical field — say about 
650 kilovolts per cm. Thus for TF = 10~ cm, or 50 kmc, Vd would have 
to be less than 26 volts. For uyL = 10 ohms, the maximum power output 
would then be less than 2 watts. At TF = 10~ cm, or 5 kmc, the maxi- 
mtmi Vd would be 140 ^'olts, which gives a power of 50 watts. At TF = 
10"^, or 500 me, the maximum voltage would be 1400 volts, and the 
minimum 100 volts. The power woiUd then be between 25 watts and 
5 kw. 

1.6 Zener Current 

Chynowcth and McKay* have shown that in sufficiently narrow junc- 
tions, where breakdown occurs at around 10 volts or less, the current is 
generated not by multiplication but by internal field emission. That is, 
electrons tunnel from the valence to conduction bands. This is known as 
Zener current. We .shall show in Section IV and Appendix E that the 
device will operate on Zener current, but much less effectively. This, 
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even more than the limited power, may limit the minimum W and hence 
the iiiaximum frequency that would' be practical. 

Expcmnents 

An experimental progi'am has been undertaken to construct and test 
a device operating liy avalanche. 

II. CUKliENT AND SPACE CHARGE 

In this section we consider the physics of multiplication, space-charge 
and carrier flow in more detail aud obtain the ecjuations that govern the 
behavior of the diode. 

2.1 MiiUiplication 

Electrons moving with velocity v generate hole-electron pairs at a 
]-ate of avn, where n is the electron concentration and a is the ionization 
I'ate. By definition a is the munlier of pairs produced on the average by 
an electron in moving unit distance. Thus a'^ is the average distance 
between ionizations. We shall be dealing with the case where each pair 
produces roughly one other pair; thus a~^ is of the order of the \\'idth 
of the narrow region near the n'^-p junction where the multiplication is 
localized. For fields below about 000 kilovolts per cm a can be con- 
sidered a function of field, a = a(J'J). For larger fields the electrons do 
not have time to get into eqnilibrinm with the field in the average dis- 
tance a"^ between colhsious. We shall be deaUng with fields where a = 
aiE). 

McKay^ has measured a as a function of E for electrons in sihcon. 
Fig. 3 is a plot of the best fit curve to McKay's data, together with the 
theoretical curve calculated by P. Wolffl A straight hne of slope 6 is also 
shown and is seen to be a good fit to the data below about 500 kilovolts 
per cm. In the formulas we shall use the general relation a ^ E'" and 
take m = G for numerical calculations. 

Assum-piions 

We shall make the simplifying assumption that a = a{E) is the same 
for both holes and electrons. Actually this is not so. However, a is so 
sensitive to field that the difference in field for the same a is small. Any 
difference in the a's will alter both the carrier and field distributions so as 
to favor the carrier with the lower a. Consequently, it is believed that 
this assumption will not give misleading re.snlts. 
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We also assume that holes and electi'ons travel with the same velocity 
in the high fields involved. This will not lead to serious error since the 
motion of one type of carrier always plays a niinor role. 

2.2 DC Case 

We take the x axis normal to the junctions with .-c = at the n"'"-p 
junction and x = W at the i-p^ junction, as shown in Fig. 4. It can be 
shown that the direct current / is related to the thermally generated re- 
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verse saturation current /, by 

^^ ^ 1 - j a(E) dx 



(8) 



wluin' tiic integral is taken over the space-charge region. See, for exam- 
ple, McKay^. We will be dealing with space-charge regions where most 
of the multiplication occurs in a narrow region, which we shall call the 
mulliplication region, near .t = 0. Thus we need take the integral in (8) 
only from x = to x ^ 3:1 , where .ri « W is the width of the narrow 
multiplication region. 

Breakdown occurs, that is, the direct current becomes infinite, when 



/ 



a (Ix = 1 . Physically this means that each hole-electron pair generated 
ill the multiplication regions will, on the average, generate one other 
pair. We now consider how i adx depends on the peak field Eo for vari- 
ous field distributions including the one shown in Fig. 4. 

(a) The simplest case is that where the field is a constant Ea in the 
narrow multiplication region. This would correspond to a short region 
of high constant field followed by an abrupt drop to a much lower field. 



1 ^o - Ec 




W 



*-~-l/(m+i) L 

Fig. 4 — Field distribution at breakdown 
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Then 

/ a dx = / a{Ei) dx - Xia{Ei). 
So if a is proportional to E'" 

where Ec the critical field for breakdown is given by a{E[) = l/a;i . 

(b) Next consider the linear field distribution shown in Fig. 4. Curve 
A is the field distribution at breakdown. We can take E = Eq — kx in 
the midtiplication region. The slope /,■ is proportional to the charge in 
the p-rcgion. We shall neglect the effect of the carriers on the space 
charge, Then /,' wUl be proportional to the fixed charge in the p-region 
and will lie constant independent of current. In practice k would be well 
above Er/W. In order for the carriers to produce a field gradient com- 
parable to this, the current would have to be of the order of (kv/^tt) 
{E,-/W), or several thousand alnperes/cm^ Actual currents will be much 
smaller than this. 

In Fig. 4 the flat section of field \vill contribute negligibly to the mul- 
tiplication, so we can replace the actual field distribution A by the 
tangent curve B. Then 



/-•»^(0 



(if) 



where Ec is given by a{Ec) ^ (m -\- l)/l and I ^ Ejk is the zero inter- 
cept of the tangent curve B at breakdown. Thus the peak field E,: at 
breakdown is equal to the breakdown field for a constant field region of 
width l/(m + 1), as shown by cur\'e C in Fig. 4. 

(c) For a linear gradient junction the field varies parabolically with 
xand 






In what follows we shall use (9), which applies to a linear field, or 
abrupt junction. This will be perfectly general if m is left arbitrary. The 
results obtained for m = G will be a good approximation for a linear 
gradient junction, since the difference between 7 and 6.5 is within the 
experimental uncertainty in m. 

As an example, consider a step junction with m = 6 and W = 10~ 
and let I = 0.717. Then the critical field E^ is given by a{Ec) = 7/1 = 
10~* cm. From Fig. 3, E^ - 350 kilovolts/cm. 
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2.3 Time Dependence 

In the dc case / a dx cannot be greater than unity. That is, we can- 
not get above breakdown. This is not necessarily so for rapidly varying 
fields. Wg now derive a differential equation for the current as a function 
of time. This will reduce to (8) in the dc case. 

Let p and n be the hole and electron densities respectively. The cor- 
responding eiu'reuts are Ip ^ qvp and /„ ^ —qvn, where q is the elec- 
tronic charge. Hole-electron pairs are being generated at a rate av(n -\- p). 
This is so large compared to the rate of thermal generation that the 
latter can be neglected. The equations of continuity then become 



t= -if' + -(« + ?), 
dt q ox 



dn 
Tt 



q dx 



av{n -\- p), 



(10) 



(11) 



w^here a = a{E). The three variables, p, n and E are determined by 
(10) (11) and by Poisson's equation and the boundary conditions at the 
junctions. An exact solution is impractical. Instead, we make an ap- 
proximation based on the fact that most of the generation is confined to 
the narrow multiplication region near x = 0. This is shown in Fig. 5. 
The multiplication region extends from x = to a; ^ Xi . The width .Ti 
is a small fraction of the width 11' of the space-charge region. 

2.4 Assumptions 

We shall assume that (a) all of the nudtiplication occurs in the rela- 
tivelj' narrow multiplication region and (b) the total current I{x, t) = 
Ip{x, t) -\- I,i{x, /) in the multiplication region is a function /{x, t) — /o(0 

MULTIPLICATION 
REGION 




Fig. 5 — Boundary conditions on the multiplication region. 
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of time only. This, of course, will not be tme of Ip and /„ individually in 
either the ac or dc cases. Assumption (a) will be valid if the net acceptor 
concentration in the p-region is sufficient to give a sharp drop in field 
and hence to localize the multiplication. The second assumption will be 
a good approximation provided the current does not vary by a large per- 
centage in the time n = Xi/v required for carriers to cross the multipli- 
cation region. This will not be so at large amphtudes. However, the 
errors for rising and falling currents tend to cancel, so the equation is 
right on the average. As we shall see, at large amplitudes only the 
average is involved. 

2.5 Differential Equation for Io{t) 

Adding (10) and (11), using the assumption /p + /„ ^ Io{t^ and in- 
tegrating from X = to 3; = 3;i gives 

n^= -ih- InV + 27o r a dx, (12) 

at Jo 

where n = xi/v is the transit time across the multiplication region. The 
boundary conditions are shown in Fig. 5. The hole current at .-c = 
consists entirely of the reverse saturation current /ps of holes thermally 
generated in the n"*"-region; these have moved to the n'^'p-j unction by 
diffusion. Thus at a; = 0, 7^ — /„ = 27p — 7o = 2Ips — h ■ A.t x = xi 
the electron current consists of the reverse saturation current 7„a of elec- 
trons thermally generated both in the space-charge region and in the 
p"^-region, so 7p — 7„ = —2I„, -\- h . With these boundary conditions, 
(12) becomes 

In the dc case 7o is the direct current 7, so this reduces to (8). 

The condition for breakdown is / a rfa; = 1. If a field that satisfies 

Jo 

this is suddenly applied, 7o wiU increase linearly at a rate of 21 Jn and 
become infinite. If a larger field is apphed, 7o will approach infinity ex- 
ponentially. For a smaller field h will approach a finite value. 

The integral I adx over the multiplication region depends very little 

on carrier space-charge. Hence it will be negligibly affected by the small 
differences in carrier distribution between the ac and dc cases. We can 
therefore combine (9) and (13) to obtain a differential equation 

^|hi7„^ {Eo/Er"-' - I + r (14) 

2 at 10 
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relating the current h = h{t) in the multiphcation region to the peak 
field £"0 = Eo{l). 

In most practical cases the current h will be so large compared to /. 
that the effect of I» can be neglected. The correction due to /, and the de- 
tiiiled formulas for evaluating the effect are given in Appendix D. 

2.G Example 

At low enough amplitudes of oscillation we can expand Eq/Ec in powers 
of Eo/Ec — 1 and retain only the linear term. Then, neglecting /, , equa- 
tion (14) becomes 

Uah^ '^"'+'H §^-i). (15) 

at Tl \tte / 



If /o iH to be periodic, then Eo must be periodic and the dc bias must 
be such that the average Eo is Ec . Suppose we apply a periodic voltage 
with the proper bias so that Eq = Ec + E^ sin ut, where a: is the opti- 
mum frequency -k/t. Then 

/o(0) TV TiEe\ t/ 

Fig. 6 shows the field and Inh as functions of time. Suppose r/n ^ 
W/xi = 20 and m ^ 6. Then even if the amplitude Ea of variation of 
En is as small as 1.8 per cent of Ec , h will vary by a factor of 10 over a 
cycle. Thus we can have small signals in field and voltage liut large sig- 
nals in tnirrent. We shall call this the intermediate range of amplitudes. 

From (Hi) the maximum value of h is seen to occur in the middle of the 
cycle, where I ^ r. Thus, if la varies by a large factor, the current is 
generated mainly in a pulse in the middle of the ac voltage cycle, as 
shown in Fig. 6. 

Actually the space charge of the current will affect the Efi{t) cur\'e, 
wliich will not be exactly sinusoidal for a sinusoidal voltage. However, 
this does not affect the conclusion that a small field and voltage signal 
can give a large current signal and that the current /« approaches a pulse 
as the ac amplitude increases. 

2.7 Carrier Space-Charge . 

We have now dealt with the multiplication region, and obtained an 
equation relating /„ and E„ . It remains to consider the rest of the space- 
charge legion, where current generation is negligible. We shall also 
neglcci Ihe re\erse saturation current /,„ of electrons, which will be 
negligible compared to the total current I{x, t). From Poisson's equation 
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we shall derive another relation between h and Eq . This will involve 
also the voltage V(t) and together with (1-4) will uniquely determine 
hit) and Eo{t) for any V(t). 

2.8 Physical Picture 

The physical picture is shown in Fig. 7(a). The width of the narrow 
multiplication region is small compared to the total width 17 of the 
space-charge region. Therefore, in treating the current, space-charge 
and field distributions throughout the whole .space-charge region, we 
shall assume that the multiplication region has zero width so that all 
current is generated at x ^ 0. Then Io(t) is a current of holes flowing 
out of X = 0, and the only carriers in the space-charge region are the 
holes. 

The current I(x, t) at any point x and time t is 

I{x, t) = 7(0, t - x/v) - 7o(( - x/v). (17) 

Thus the entire hole and current distributions are given in terms of 
/o(/). 




1 




Fig. 6 — Case of a sharp current pulse. 
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Next consider the efteet of the holes on the sp:\ce charge and the field. 
The field distrilnition PJp^x) at punch-through is shown hy the dashed 
curve in Fig. 7{b). If no current flowed the field distribution for any 
voltage above the punch-through voltage would he simply A'„(a) plus a 
constant determined by the voltage. The fixed negative charge in the 
p-material gives a drop in field across the space-charge, or depletion, 
layer. The positive space charge of the hole.s opposes that of the ac- 
ceptor.^ and hence reduces the drop in field. So, for a given voltage, £'0 
will decrease as the current flowing in the space-charge region increases 
and flattens the field distribution. 

2.9 External Current 

The holes traversing the .space-charge region give rise to a current 
/, = /,(/) in the external circuit. /, is e(iual to the average current flow- 
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Fig. 7 — r'ljrreiit Ui) and field (h) in the s|jace-('h;irge region. 
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iug ill the space-charge region. (A rigorous and simple proof of this for 
the present case of plane parallel geometry is given in Appendix A.) Thus 

■ hit) = ^ j I{z, t) dx 

^ ^^ j\.{t - x/i^ dx (18) 

= - f' Ut') dt'. 

T Jt-T 

In other words, the current le in the external circuit is the total charge in 
the space-charge region divided by the transit time t. Fig. 6 shows le for 
the case discussed at the end of Section III, where h was a sharp pulse 
in the middle of the cycle. From (18) it follows that the average value of 
Ic(t) is equal to the average of Io{t). 

In addition to h , which arises from carriers moving through the space- 
charge region, there is a capacitati^'e current 

dV 
A-^l^ (19) 

flowing in the external circuit. When the voltage V{t) across the diode is 
above the punch-through voltage, Vp , the capacity is a constant given 
by (2). As discussed earlier, Ic is the current required to charge and dis- 
charge the diode regarded as a capacitor. It furnishes the variation in 
charge at the edges of the space-charge region. 

2.10 Effect of Current on Field 

We now show how the space charge of the holes reduces the peak field 
for a given voltage. The stability of the device comes from the fact that 
current multiplication increases as Ed increases but the current carriers 
give a space charge that reduces En . 

If there were no current flowing, any increase in Y above Vp would 
simply raise the entire field di.stribution by an amount {V -~ Vp)/W. 
Fig. 8(a) is a plot of the difference E'{x, t) = E{x, i) — Ep(x) at a given 
time. The slope of the £"(x, t) curve is determined entirely by the space 
charge of the holes; the effect of the fixed charge is already included iu 
Ep{x). The holes give a charge density I(x, t)/v. Thus, Poisson's equa- 
tion is 

^ {x, f) ^ ^ I{x, i) ^ ^ hit - x/v). (20) 

dx KV KV 
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The excess of V over Vp is equal to the area under the curve in Fig. 8(a). 
This is equal to WE'{0, t) = ^V[Eo{t) — Epo] plus the sum of the areas of 
a number of horizontal strips like the one shown. The area of such a strip 
is (IF - x) dE'. So 



(IF- ;r)--rf.r. 
(1 o.i' 



Substituting (20) into this and setting t' = t - x/v gives 

no - 1\ 



Eo{t) - Epo + 



TT^ 



+ AE(t\ 



AE{t) = -~ f hiOW -t-\-t'] dt'. 

KT Jt-T 



(21) 



Here AE(t) is the effect of current on the field. The quantity AE is al- 
ways negative. 



E'(X,t} = E[X,t)-Ep[X) 




E 



Eo(t)-Epo 



(b) 



w — ♦x- 



Fig. 8 — Effect of carrier space-chiirge on field. 
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2.11 Effect of a Current Pulse 

To illustrate the physical meanmg of (21) we return to the example 
illustrated in Fig. 6. A sharp pulse of current was generated near the 
middle of the cycle. The reduction —AE in Ea due to this pulse is also 
shown in Fig. G. If the pulse is instantaneous, —AE jumps at once to its 
maximum value and then declines to zero linearly in the time r required 
for the pulse to cross the space-charge region. The physical reason for 
this is easily seen by reference to Fig. 8(b), which shows the same thing 
as Fig. 8(a) except that the carrier space-charge is concentrated at one 
point, that is, in a pulse. For a given voltage (area under the curve), the 
reduction in Eo will deci'ease from its maximum ^'ahie to zero as the pulse 
moves to the right across the space-charge region. From the effect of a 
single instantaneous pulse of current on Eo , the effect of any arbitrary 
current pulse can be found by resolving it into a series of instantaneous 
pulses and superposing the effects. An instantaneous pulse of current of 
charge 5Q give.s an instantaneous drop in field of 4x5Q/k. 

III. ANALYSIS 

We now have two equations, (14) and (21), relating the current foil), 
the field Eo(i) and the voltage V(t). Thus in principle the current can be 
found for any applied voltage. Actually the exact solution is impractical 
except in limiting cases. In this section w^e present (a) an exact sohition 
for the linear .^mall-signal case, (b) an approximate analysis for large 
amplitudes and (c) a rapidly converging iteration method that yields 
solutions of any desired accuracy. 

3.1 VuUage 

We shall assume tliat the voltage varies sinusoidally, V{t) = Vd + 
Va sin io/. This will certainly be a good assumption in the small-signal 
range, where the diode is hnear. The cavity is hnctir at all amplitudes. 
At large amplitudes of oscillation a sinusoidal voltage gi\'es a sinusoidal 
capacitative current Ic plus a conductive current Ic which approaches a 
S(|uare wave as the amplitude increases. Thus we are assuming that the 
voltage across the cavity is sinusoidal while the current contains a dis- 
tribution of higher frequencies. For a cavity Hke that shown in Fig. 2, 
this assumption may be a relatively good approximation. If the cavity 
is tuned to the fimdamental it may be almost a short circuit for the 
higher frefiuencies. 

3.2 DimcHsiiitilcss ]'arinbl('s 

It will simplify the discussion to express everything in terms of dimeu- 
sionless variables, taking as units parameters that characterize the de- 
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vice. For cxamplp, IF would be unit length, t, unit time, find Ec , unit 
field. Then the carriers would travel with unit velocity v = W/t = 1. 
As can be seen in Fig. 1, the voltage will be of the order of, but less than, 
the unit voltage WEc . It is convenient to choose the unit charge so that 
^tt/k = 1. Then unit charge produces a unit gradient of field. Since 
}) = 1, a unit charge moving in the space-charge region gives unit current. 
Hence unit hole current in the space-charge region produces imit field 
gradient, and the a\'erage current le in the space-charge region is equal 
to the total drop in field due to the carrier space-charge. The aciual ciu-- 
rent will be small compared to unit current since the curent produces 
a drop in field that is smaller than the ac field variation, which in turn is 
small compared to Ec . 

From the choice of units it follows that the diode has unit capacity 
and that optimum frequency is to ^ x. The following table sununarizes 
the units and gives typical \-alues for a silicon diode with TF = 10' cm. 
McKay's data,' plotted in Fig. 3, gi^'es A', = 350 kilovolts per cm; the 
effective width of the multiplication region at breakdown is taken to be 
10^* cm. 

Quantity Unit Example 

length ir 10-= cm 

time r lO-'o sec 

field E, 3.5 X 10* volts/cm 

voltage WEa 350 volts 

current density ^ = vCEc 3.7 X 10'' .im]>.s/(^iii= 

■iTTT 

]30wer (li^nsity £ £.= 1.3 X ]()« wLifts/r-m^ 

The unit of power is seen to be relatively independent of W, since unit 
voltage goes as W and unit current as 1/lF. (A given current causes a 
greater drop in field across a wider space-charge region.) 

3.3 Governing Equations 

Since the peak field Eo{t) varies around Ec it is convenient to define a 
dimensionless field 

^(0 = m-i^, (22) 

Then, for a sinusoidal voltage variation, equation (21) becomes in dimen- 
sionless terms 

E{t) ^ Eh + T'. sino^f + AE, 

r' (23) 

AE = - h{l')[l - t + t']dt'. 

Je-i 
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Here A^" is the effect of current on E and Eb the effect of the dc voltage 
Vd • Let Epu be the peak field at punch- through. Then, at punch- through, 
E is (EpQ — Ec)/Ec . If no current flowed, any increase in dc voltage above 
the punch-through voltage Vp would simply raise the whole field dis- 
tribution by (Vd — yp)/W (see Fig. 7). Thus the dimensionless parame- 
ter Eb is given in terms of dimensional quantities by 

The value of Eb will be very small compared to unity. 

In equation (14) we neglect /g and expand the right-hand side in terms 
of powers oi E = E{i). It will be sufficient to stop the expansion at ^. 
Then 



^ In 7o = ^i?^L±i) [E+'^g'). (25) 

at Ti 



The current through the diode is the capacitative current 

L = C dV/dt - oiVa cos (^t (26) 

plus the conductive current /« , where, from (18), 

IM = f hit') dt'. (27) 

Average, or dc, Values 

We shall let ( ) denote time averages and define Id as the direct cur- 
rent Id = (le + Ic) ^ (le) = (lo) where the last step follows from (27). 
Averaging (23) gives 

(E) = E,-^ (28) 

or (A£?) = —Id/2. Since h{t) is periodic, we have, from (25) 

{E) -f I {E') ^ 0. (29) 

Thus, in the small-signal limit the average field {E) vanishes, and Id = 
2Eb , where Eb is given by (24) . 

3.4 Linear, Small-Signal Case 

When all quantities vary by small fractions of their average values, 
then the equations separate into a dc part and a linear ac part which is 
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easily solved. We now derive the impedance as a function of the direct 
current Id for the optimum frequency ta = ir. The ac voltage is Va sin -wt 
or, in complex form, "Tac'""'. We write h = Id + /oae'""' where 7o„ is com- 
plex. Similar expressions can be written for le and E. Equations (23), 
(25) and (27) give the ac relations 

2(m + 1) 



tV/oo = 



ri 



lir 



ha - ^ /., 



Id^a , 



- Ir 



(30) 



2^.1- 






(31) 



iV(T'„ - 7i'„) ^ /o 
From equations (30), Va = Zha where 

2 L 2(m -{- l)/jj 

The current through the diode has two parts, Ic and h , where 7^ is 
a pure capaeitative current. Thus the equivalent circuit consists of a 
unit capacity in parallel with an element of impedance Z, where 7^ goes 
through Z. The impedance Z consists of a fixed reactance and a resistance 
that varies with the dc bias. A simpler equi-\'-alent circuit emerges from 
considering the admittance Y of the diode. Since the capacity and the 
impedance Z are in parallel^ )' ^ itt -\- 1/Z. Here the Q of the diode is 
the riitio of the imaginary part Yi of Y to the real part Yr . From (31) 



^ 2 



1 - 



IT n 



2(?n + l)7j' 



Yr = 



tQ 



(32) 



F, = 



Thus Q varies hnearly with Ti/Id and is negative for 7^ less than 
[T^ri/2(??i + 1)]. When Id is equal to this critical value the diode becomes 
an open circuit for this frequency. This means that none of the alter- 
nating current generated in the midtiphcation region flows out of the 
diode. Rather it flows into the edges of the space-charge region and pro- 
vides the current that charges and discharges the diode regarded as a 
capacitor. In other words, the unit acts as a capacitor generating its own 
charge internally. Hence voltage can vary with no external alternating 
current. 

3.5 Equivalent Circuit 

E(juations (32) describe a simi}lc equivalent circuit consisting of a 
fixed unit capacity in series with a conductance irQ, whei'e Q depends on 
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the dc bias. The equivalent circuit suggests the followmg practical con- 
clusion; If the cavity is designed to act like an inductance in series with 
a variable load resistance, then the load resistance can be made equal to 
the negative resistance of the diode, so that the two resistances cancel, 
and the equivalent circuit consists of the fixed capacity of the diode and 
the inductance of the cavity. Hence at small amplitudes the resonant 
frequency will be independent of the dc bias. 

3.6 Sharp Pulse Approximation 

As the oscillation builds up, the behavior rapidly ceases to be linear, 
and it is impractical to solve the equations exactly. However, as we have 
seen, when the amplitude increases the current Io{t) approaches a sharp 
pulse and Ie(t) approaches a square wave as shown in Fig. 6 for the opti- 
mum frequency w = tt. The a.\'eragc current Id is half the maximum. In 
the limit of a sharp pulse the problem again becomes simple. We now 
derive some approximate relations for this case, and show how the 
oscillation can be stabilized. 

The power delivered Pr is 



p,^ ''■ 



~ I Icii) sin irt dt. 



We call sub.stitute (27) for IJ4,) into this and reduce the double integral 
to a single integral l)y integrating by parts. The result is 

Pr = - '—{h{t) oosttO (33) 

TT 

where again the brackets denote the average over a cycle. 

Thus, if the current /o is generated in a pulse near the middle of the 
cycle, where cos irt is negative, Py will be negative and power will be de- 
livered to the ac signal. 

We define the Q of the diode as 2t times the ratio of the energy stored 
in the capacity to the energy lost in a cycle. The stored energy for unit 
capacity is F//^- The energy lost is the negative of the power delivered, 
Pt , times the period 2T/aj = 2. So 

1 ^ £ (7o(0 cos irt) , , 

Let the pulse of current occur at a time ti . Then in the hmit of an in- 
stantaneous pulse (34) becomes 

7^ ^ —jT cos xii . (35) 

Q ir~ Va 
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The phase relations will be ideal for /j =- 1; that, if^, when the pulse 
occur.s in the middle of the cj'cle. In this rase —Q iiiei'eases with T„ for 
constant direct current, so the oscillation is stable.* 

We now consider how h depends on T'„ and the dc bias. The current 
pulse becomes sharper as Va increases and n decreases, as can be seen 
from (IG). In the limit of vanishing n the pulse becomes instantaneous 
and the problem can be solved exactly. This is done in Appendix B. 
Here we give a simple physical argument that will be a good approxi- 
mation .so long as the duration of the pulse is small compared to a period, 
as will be the case in the range of practical operation. 

Kig. 9 shows E{t) = Kb + T'„ sin iri + A7i' for the case where the 
current pulse occurs at an arbitrary time h ■ As illustrated in Fig. (i, a 



E(tl 




Fig, 9 — Variation of peak field witli time for sharp current pulse. 

current pulse causes ^E to drop abmptly and then rise hnearly to zero 
in a transit time. Let Ei be the abrupt drop in E. We have already seen 
that the average i^E is -Id/2. Since A-E is a triangular wave, lasting 
half a cycle, it follows that the maximum is four times the average; so 
E^ =^ 2/j . 

Let us consider in some detail the relatively short uiterval during 
which current is being generated. The holes can move only a short dis- 
tance durhig this time. Hence, as seen from Fig. 8, the field will drop 
by an amount roughly proportional to the amount of charge generated. 
It will, therefore, have dropped by E/./^ when half the charge has been 
generated. The pulse will be roughly .symmetrical and will have I'eached 
its peak when half the charge has been generated. The current is a maxi- 
mum at E(t) ^ since it builds up for positive E and decreases for 
negative E. Thus, during the first half of the pulse E drops to zero, the 



* I !im indel>ted to J. L. Moll for pointing out the advantage of applying the 
bia.s with a eoii.'^tant direct current, generator. 
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drop being E&/2. During the second half of the pulse the held continues 
to (h'op to —E^I'l. This is shown in Fig. 9. Before the pulse oecurs the 
held is Fj{() = Eb + sin -kI. Hence in the limit of an instantaneous pulse 

E^, + F„ sin -Kk = Ei/2. Since Es - 2Li 

Id = Ef,+ F„ sin irk . (36) 

Ehminating {E) from equations (28) and (29) gives another relation 
between 7^ and Eb : 



/., = 2R, + (f)<B'). 



Since Eb , Va and — A^ ^ 21 d are ah small compared to unity, the only 
contribution to {E') that cannot be neglected in comparison with Id 
and Eb is {Vl sin' irt) = Va/2. So we have 

/, = 2Eb + ("^ V:. (37) 



Ehminating Ei between (3G) and (37) gives 

sinxi, ^ A+|k. (38) 

When the right-hand side is larger than unity the current cannot be 
an instantaneous pulse. What happens then is that the current varies 
almost sinusoidally and produces a space charge which keeps E small 
at all times. In this (;ase the phase relation between the voltage and cur- 
rent makes Q positive. 

In practice, we begin with a small enough bias current so that Id/2Va 
has become small before the oscillation has built up to the range where 
the sharp-pulse approximation is valid. The effective Q of the cavity will 
be chosen so that the oscillation will be stabilized before (m/4) Va be- 
comes comparable to imity. For example, suppose we take Id ^ Fa/2 
as discussed in Section I. Then —cos tt/i will have dropped only to 0.85, 
even for as large an amplitude as Va = 0.2. For Va = 0.1, cos irti would 
be -0.92. 

3.7 Constant dc Voltage Bias 

If a constant dc voltage bias Vd is apphed then Eb is constant. As the 
oscillation builds up Id increases in accord with (37). 
In terms of Eb and 7„ , equations (35) and (38) become 
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and 



, Eb . m j^ 
sm Trti = ~ -\- — Va . 



(40) 



The solid curves in Fig. 10 are plots of —l/Q vs Va for two values of 
Et . In the range shown, cos irti is approximately —1. The curves, there- 
fore, have minima at Va = 2\/Eb/m- As Va increases, —l/Q will even- 
tually reach a maximum and begin to decrease because the term (?n/2)Fo 
in (40) will become important, and the pulse will occur too soon in the 
cycle. This range lies beyond that shown in Fig. 10. 

At the amplitudes shott^l in Fig. 10, and especially in the stable range, 
the main error in the above approximations comes from the fact that 
the pulse is not sharp. We now turn to a method of obtaining more ac- 
curate results at low amplitudes. 

3.8 Iteration Method 

Equation (25) gives /o(0 when E{t) is known and (23) gives E{i) 
if /o(0 is known. Thus we can guess at E{t), find the corresponding 
h{t) from (25) and use it to determine a new E{t) from (23), and so on. 
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Fig. 11 — Converging solutions for /o(0- 



For the initial E(t) we take the instantaneous pulse solution. The pro- 
cedure converges rapidly and results in solutions of any desired ac- 
curacy. In finding /«(/) from E{t) in (25) a constant of integration, 
/o(0), is involved. What is found directly from E{t) is h{()/hiO); /,i(0) 
is so chosen that the next E{t) will satisfy the condition (29) that the 
following foil) be periodic. The procedure is discussed in more detail in 
Appendix C. 

Fig. 11 shows plots of hit) for n - 0.05, E„ - 0.004 and T,. = 0.0;i. 
The various itei'ations ;i.i-e numbered. The values of —\/Q in successive 
iterations were (to four places): O.UKi, 0.1330, 0.1321 and 0.1321. Fig. 
12 shows — A/i' vs / for n ^ 0.05, Ei, -= 0.004 and various values of 
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Va . The coiTcsiJoiiding /o(/.) curves are shown in Fig. 13, and Jt(t) 
curves are shown in Fig. 14. 

Fig. 1 is a plot of —\/Q vs Va for various values of E^ and n . The 
points were obtained from the iteration procedure and the solid curves 
from the sharp-pulse approximation. As expected, this approximation 
improves as Va , Eb and l/n increase. With Eb ^ 0.004 there is a mini- 
minii for ti = 0.05 but none for ti ^ 0.1. In the former case stable os- 
cillations would he possible at ampUtudes below about 0.035. 

3.9 Inkrmcdiaie A mplitudes 

We now have solutions for the hnear small-signal range and for the 
large amplitude range, where the sharp-pulse approximation is good. 
However, there remains an intermediate range where the current varia- 
tion is too large for the small-signal analysis to apply but not large enough 
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■for the sharp-pulse approximation to be vahd or even practical as a start 
for the iteration procediire. This range can be dealt with by an approxi- 
mation that becomes more accurate as — Q increases and can be within 
10 per cent even for —Q as low as 7. The results show how /o(0 changes 
from a cosine wave to an increasingly sharp pulse as Va increases. 

If Q~^ is small compared to unity then, in the small-signal v,ase, ~Q = 
irVJIea , where /f„ is the amphtude of variation of h . The same is 
true at large amplitudes if I^it) is analyzed into a Fourier series and 
lea taken as the amplitude of the fundamental. In the linear small- 
signal range, I^ is proportional to Va so — Q is constant. If I4 is kept 
constant as the amplitude builds up, then 7™ increases less rapidly than 
Va ■ This can be seen from the fact that V„ can increase without limit 
while Ic{f) approaches a square wave of amplitude Id , for which /,« = 
(4/T)/d . Thus Q approaclics - {jr74) (F„//d) as given by (35) , for ^ = 1 . 
The phase shift, vt: , between voltage and le is cot"' Q. We have seen 
that Ic is always 180° out of phase with the peak field. Therefore, if —Q 
is large enough so that /« and voltage are approximately 180° out of 
phase, the peak field is in phase with the voltage and we can use equa- 
tion (16) with Ea = Va ■ The current Io{t) is then 

hit) = iMe-^'-""''', 

2(m+ 1)F„ (41) 



X = 



TTTi 



This is seen to approach a pulse of increasing sharpness as !'„ increases, 
and to reduce to the small-signal results as Va approaches zero. 

The Q can be found from (35) and (41). From (35) 1/Q = (4/t ) 
(VF-.)/(i-), where 

\h{i)} ax Jq 

The function /(.i;) is (he ratio of the first to the zero order Bessels func- 
tions of pure imaginary argiunent. As x decreases, f(x) approaches x/2, 
which means that Q is constant. As x increases, f(x) bends over and 
approaches unity asymptotically. Values are /(I) = 0.45, /(2) = 0.70, 
/(4) - 0.87, /(8) = 0.94 and /(16) = 0.97. When h is constant, 
Qo/Q ^ 2f{x)/x, where Qn is the small-signal Q. Thus the curve of — l/Q 
vs Va starts out flat and then decreases and approaches the form \/Va 
as Va increa.'ies. This is illustrated by curves 1 and 2 in Fig. 15. 

If the — 1/Q vs Va curves for constant Id are calculated from the sharp- 
pulse approximation [equations (35) and (38)], they have the form of 
the solid curves 3 and 4 in Fig. 15. As T'„ decreases, the calculated — 1/Q 



432 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1958 

goes through a maximum. For Id below about 4/m, the maximum would 
occur at about V^ = Idf\/2. However, as Fa decreases, the sharp-pulse 
approximation will break down. If h is small enough ko that the small 
signal Q is negative, then the sharp-pulse approximation must break 
down before the maximum is reached, since, as we have seen, ~1/Q 
increases monotonically as y„ decreases. Ro the cin-ves have the form of 
curves 1 and 2 in Fig, 15. However, if the small-signal Q is negative, 
then the curves have the form of curves 3 and 4 and the sharp-pulse 
approximation breaks down in the lange of positive Q (shown dotted in 
Fig. 15). In this range the current varies roughly sinusoidally and pro- 
duces a space charge that keeps the field variation small. 

IV. OPEHATION 

In this section we consider in more detail some of the practical ques- 
tions about the design and operation of the diode. In particular, we dis- 
cuss the stability for both constant current and constant voltage bias, 
the limitations on both dc bias and ac amplitude, the eft'ects of heating 
and finally the frequency dependence of the effective admittance. 

4.1 Stahilily at Constant Direct Current 

Fig. 15 shows the form of the variation of -1/^ with voltage ampli- 
tude Va when the bias is applied with a constant (lii'e<'t cuiTent generator. 
The oscillation will be stable at any point where -l/Q is decreasing 
with Va and the (^'s of cavity and diode are equal in magnitude. The 
horizontal dashed line in Fig. 15 represents the effective \/Q of the 
cavity. Thus the oscillation can be stabiHzed at the points A, B and C. 
The curves are nimibered in order of increasing direct current Id . If 
the direct current is turned on slowly compared to the response time 
Qt of the diode, the oscillation will begin when h is slightly above the 
value for curve 1. Thereafter the oscillation will build up and V^ will 
increase as Id increases. By slowly varying Id we can establish a stable 
oscillation at any desired amphtude. 

By raising h rapidly enough we could get onto curves 3 or 4 while 
the amplitude was so low as to be in the imstalile range. The oscillation 
would then start but quickly die out. 

Because of the nonhnearity of the diode, the dc voltage will vary as 
Va increases at constant h ■ Ehminating Eb between (24) and (37) gives 

/,; = 2(F. - Y, + E„n - 1) + ^ \\\ (43) 
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where A'po i« the peak field at punch through (taking E^ as unit field) 
and hke V p is a couwtant of the diode. 



4.2 Stability at Constant dc Voltage 

When the bias is applied by a constant voltage, Et will be constant 
as seen from equation (24). Curves for several values of Eb and 
Ti were shown in Fig. 10. The small-signal Q is given by (32) with 
Jd = 2Eb . As Eb is rai-sed, the small-signal Q will become negative 
and — 1/Q will rise. If Qc , the effective Q of the cavity, is snuUl enough 
the oscillation will initially be stable. For example, in Fig. 10, if 
Ti = 0.05 and l/Q^ = 0.10, the oscillation will begin when _£"(, is raised 
to about 0.004. As Eb is further increased, the amplitude !'„ will 
also increase. However, as seen from the figure, the oscillation will 
not remain stable as Eb increases unless l/Qc also hicreases. For 
example, at Eb = 0.008 the minimum —l/Q is almost 0.18. We have 
seen that at lai'ge enough amplitudes —l/Q reaches a maxinmm and 
decreases to zero. Thus, if E^ is increa.sed with constant Q^ , the ampli- 
tude will sud<lenly jump from the range shown in Fig. 10 to a much 
larger value. In the range where the sharp-pulse approximation is good, 
— l/Q is given l)y (37) and (40). For example, when Va is large com- 
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pared to {^/m)Eb , we have, approximately, 

1 4m ^. . „ 

^ " (44) 

, m -, . 

sin Ttl = — Va- 

So there is a maximum m - \jQ at W = \/2/*». The maximum will 
he at lower ampUtudes for larger Eh . 

It is seen from (39) and (40) that the oscillation will always be stable 
when Va is between t\\-o and four times E,./7H, pro\'idod Eb is large 
enough that the sharp-pulse approximation is good in that range. Thus, 
with constant dc voltage, the oscillation can be st!\bilized at any ampH- 
tude, but to reach the operting point it may be necessary to vary the Q 
of the cavity as the bias voltage is increased. 

In the remainder of this section we shall consider several effects, such 
as heating, wliich may cause the power and - 1/Q to tlecrease with in- 
creasing amplitude. These may limit the maximum power output (es- 
pecially at low frequencies) and may also be used to .stabilize the oscil- 
lation in the case of constant dc voltage. 

4.3 Efficiency 

We define the efficiency, e, as the ratio of the ac power, Pr , to the power 
Pd delivered by the dc voltage or current source; therefore, 

Pd = Vah . (45) 

The power Ph that goes into heat is the difference Ph = Pa — Pt . 
Therefore, 

P. = ^ = -^. (46) 

e 1 — e 

Later in this section we consider the temperature rise caused by Ph ■ 
At the optimum frequency, w ^ tt, 

-El = - ?Z." <^o(0 cos Tct) ,^7^ 

At small and intermediate amplitudes this can be evaluated by (41) 
and (42), which hold when -Q is well above unity. In the sharp-pulse 
approximation 

Thus, for Va ^ Vd/2 and h = 1, the efficiency would be over 30 per cent. 



e^ -^^COSirU. (48) 
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4.4 Limitations on Bias and Amplitude 

For a desired operating amplitude \vc want Yd to be as wniiill as pos- 
sible to maximize the efficiency. However, as discussed iu Section I, 
the minimum Yd is limited by the requirement that the field in the in- 
trinsic region must not become negative during the negative half of the 
voltage cycle, and the maximum Ya is hmited by the necessity of keeping 
the multiplication localized. We now consider these requirements in 
more detail. 

In the analysis, we have assumed that the field in the intrinsic region 
is always high enough so that the carrier velocity remains constant, 
independent of field. However this wdll not be so at large enough ampli- 
tudes. The field may. Id fact, momentarily become negative. This would 
reduce the negative resistance and eventually destroy it. The effect of 
field on velocity gives both an upper limit on the allowable amplitude 
for a given bias and a method of stabihzing the oscillation at any desired 
amplitude. 

The minimum field will occur in the intrinsic region and at the trail- 
ing edge of the pulse of holes advancing to the right, as illustrated in 
Fig. 8. The minimum field wall fall below the constant-velocity range 
only at liigh enough amplitudes so that we can use the instantaneous 
pulse approximation. The holes will then be closely bunched as in Fig. 
8(b). The pulse causes a drop in field equal to the current I^ , which is 
2/,i during the half cycle that the pulse is moving. The field i?,- immedi- 
ately behind the pulse is, therefore, 

E,{t) = Yd - Vp + Y. sin irt - 2/rf(l - t + ti). (49) 

This has a minimum at a time t^ where 

2 Id 

cos Xi2 = — - ^ - ,^„. 

T V„ (50) 

TT < 0)^2 < 3Tr/2. 

Equation (38) gives h , and (43) gives U in terms of Yd — Y^ and Fa . 
Thus for a given diode the minimum Ei = Ei{t2) is determined by the 
voltage bias and ampfitude. 

In general, the a\'erage field will be well above the range where ve- 
locity depends on field. Hence w^e can take -£.(^2) = as the critical 
condition where the drop in velocity in the negative half of the cycle 
begins to reduce the power appreciably. As Ya increases, the minimum 
Ei decreases. If the direct current is held constant, so that the oscilla- 
tion is stable, then we will want Yd — Vp to be large enough to liuild up 
Ya to the desired amplitude without reducing E, below zero and thereby 
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losing power. However, Vd .should be made as low as posKihle io maxi- 
mize the efficiency. Thus Eiiij) = is an optimum (condition for opera- 
tion at constant direct current. On the other hand, if we wish to operate 
with constant dc \'olt^ige and in the range whi<-h would be unstjible for 
constant velocity, then Vd should be low^ enough so that the velocity 
variation will come in at the desired amplitude and stabilize the oscilla- 
tion; so the condition EiQ^) = would not only be optimum but neces- 
sary. 

We have seen that £",(^2) depends on Vd — Vp , Va and the peak field 
EpQ at pun eh- through. Thus for a given diode Ei{U) = gives a relation 
between Vd — Vp and Va ■ The quantity Va and the corresponding 
Vd — Vp could then be chosen to maximize the power. However, an 
upper limit on the voltages is set by the necessity of localizing the multi- 
plication. The field throughout the mtrinsic region must be well below 
the breakdown field, especially when the current is flowing, r'rom Ing. 
8(b) the field at the leading edge of the current pulse is Va — Vp + Va 
sin tri + 2Id{ii — t). This determines how many hole-electron pairs the 
pulse of holes will produce ui moving across the space-charge region. 
The holes produced join the pulse and add to the power. However, the 
electrons moving in the opposite direction will disturb the phase rela- 
tions. 

In the arginnents above Vd occurs only in the combination Vd — Vp . 
It is therefore desirable to make Vp as small as possible so that Vd is 
small and e large. 

Both the voltage Vp and the peak field Epo at punch-through are de- 
termined by the impurity distribution in the p-region. The amplitude 
at which the diode is designed to oscillate will fix the choice of Epo . 
The larger the desired amplitude, the smaller Epo . However, if Epn, is too 
small the multiplication will not be localized. 

4.5 Exam-pie 

We may illustrate the above discussion by applying it to an actual 
design. In the di.scu.ssion at the end of Section I we took Vd = 2T'a = -fA/ 
as a reasonable operating condition. This would give an efficiency of 
about 30 per cent. To operate at an amplitude of Va = 0.2 would then 
require that Va — Vp be no less than 0.275. For W = 10"^ and a unit 
voltage of 350 volts this would mean Va — 70 volts, Vd = 140 volts and, 
for optimmn operation, Vp — 44 volts. These parameters are reasonable. 
In fa(^t, both Vd and Vp could probably be slightly lower so that the ef- 
ficiency would be higher. 
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4.G Heating 

The power Ph that goes mto heat will cause the temperature of the 
(Bode to i-ise above that of the surroundings. We shall assume that the 
surface of the eavity is kept at constant temperature. If the radius li 
of the diode is small compared to the radius of the center post of the 
cavity, as in Fig. 2, then heat can flow away from the diode in almost 
nil directions. The diode can be made thin enough and mounted suffi- 
ciently close to the metal so that the temperature drop in the silicon is 
small compared to that in the metal. Let AT be the difference between 
the diode temperature and the temperature of the surface of the cavity. 
Then AT is i-elated to I\ and R by the formula for spreading resistance 

Ar = ^, (51) 

where 7C is the thermal conductivity of the metal. For copper, 4X = 10.7 
watts per cm per °C. 

The temperature will rise to its equilibrium value in a time of about 
R'/D, where D is the coefficient of thermal diffusion, which isaboutunity 
for copper. 

We may now apply these results to the examples discussed at the end 
of Section 1 assmning an efficiency of I so that Pk = '2Pr . The 50 watts 
of power output at 5,000 megacycles and R ^ 0.03 cm would produce a 
temperature rise of about G0°C. At 500 megacycles and R = 0.8, how- 
e\-er, the 5 kw maximum power would raise the temperature l>y about 
600°C. Thus at low frequencies the maximum power output in continu- 
ous operation would be limited by heating rather than by how small 
the impedance w/v of the cavity can be made. However, the time con- 
stant of the temperature rise for R = 0.3 cm would be almost a tenth 
of a second, so the temperature rise in pulse operation would not be 
serious. 

4.7 I'lffcct of Temperature on Critical Field 

JMcKay' has found that the critical field increases with temperature. 
For critical fields between 250 and 500 kilovolts per cm in silicon, a 
change in tcmperatvue changes the critical field by 0.05 per cent per °C. 
Increasing the critical fieltl effectively decreases the dc voltage and cur- 
rent. Thus when the diode begins to get hot it effectively reduces its 
bias. The heating will therefore stabihze the oscillation. 
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4.8 Effect of Reverse Saturation Current 

In eilicon the reverse saturation current will become important only 
at large amplitudes. In the small-signal range it is easily shown that 
Is can be neglected if it is small compared to 7r[7rri/2(m + 1)]^ which 
will be of the order of 10"'', or a few amperes per cm for W = 10~^ cm. 

Even at large amplitudes /« will be small compared to the average cur- 
rent and will have a negligible effect on the space charge. However, as 
Va increases, the ratio of the maximum and minimum values of h (which 
varies exponentially with Va) will increase much faster than the maxi- 
mum /o . Thus the minimum 7n becomes very small at high amplitudes. 
It cannot, however, fall below /, . 

If the effect of /, is included, the equation (25) for current generation 
becomes 



^'|i„/.M^ + i)(i.' + fi.^) + ^. 



(52) 



As Va increases and the minimum 7o decreases, the term Is/In becomes 
important near the current minimum, and prevents the current from 
becoming too low. Thus /, will be important in the equation only near 
the current minimum. However, by increasing the minimum h , L 
also increases subsequent values of 7o and hence increases the associated 
carrier space-charge. The carrier space-charge wUl therefore shut the 
current off (by reducing the field) earlier in the cycle. This will reduce 
the delay between voltage and current and so reduce the power. 

The amplitude at which L becomes important can be roughly esti- 
mated as follows: From (41) the ratio of maximum and minimum values 
of /o is related to Va by 

4(m -I- 1) Inxin 

A formula for /„,,„ is derived from the instantaneous-pulse approxi- 
mation in Appendix B [equation (63)] 11. In practical cases, where la 
will be of the order of 10~', I nmn- is seen to be of the order of unity. At 
room temperature the reverse saturation current in silicon can be made 
less than a microampere per cm , which is a dimensionless current of 
around 10~^ in the kmc range. Thus the minimum /o will be greater 
than la if Va is no larger than about 2.6 n . To operate at an amp- 
litude Va = 0.2 we therefore want n to be no smaller than 0.07. 

The effect of 7. is treated in detail in Appendix D and more general 
relations are found to replace (37) and the equations derived from it. 



A PROPOSED HIGH-FREQUENCY, NEGATIVE -RESISTANCE DIODE 439 

4.9 Frequency 

So far we have considered only the optimum frequency oj — tt/t = ir. 
The small-signal analysis is easily carried out for arbitrary frequency. 
The admittance Y of the diode is found to be 

4a)(/[(l — coso)) + i(2y — sin w)] 

2;/- — 2?y sin oj + 2(1 — cos oi) 

(53) 



r(co) = 



'^^2 



1 - 



2(m -j- l)h 



} 



This reduces to (32) for oi ^ r. When ~y is well above unity, it can 
be held constant by varying Id in proportion to u. Then — 1/Q will vary 
only with the phase factor 1 — cos w. Thus by varying 7^ the device 
can be tuned mechanically over a frequency range extending from ^ 
to -f of the optimum frequency. 

At large amplitudes the current no longer varies sinusoidaUy. However, 
as mentioned in Section III, the cavity may be almost a short circuit 
for frequencies higher than the fundamental. Hence, to evaluate the diode 
as an element in the oscillator, we can analyze the current /„ + /c into a 
Fourier series and retain only the fundamental. The relation of this to 
the ac voltage then defines the conductance, G, and capacity, C, of 
the diode. The results are 

2 
G = ^nr- [{hit) shioit) sin w + (lait) coswt){l — cosco)], 

""^^ (54) 

2 

C = 1 + —^7- [{Io(i) sina:i){cosw — 1) -|- (/o(0 cos wi) sin w] 

where again the brackets denote time averages. When current and volt- 
age are approximately 180° out of phase, we can use (42) to evaluate 
the time averages in (54). 

In the sharp-pulse approximation the averages become la sin o)ti and 
Id cos w/i respectively. Equation (37), which is independent of frequency, 
remains valid. Since the drop Es in field is (47r/a;) times the average 
{—AE) = /d/2, the general form of (30) becomes 

-Id = Kb + Va sin coti . (55) 

(a 

Abo\'e the small-signal range the admittance Y = G -\- luC of the 
diode depends on voltage amplitude as well as on frequency and direct 
current. For stable oscillation the admittance Y = Y(o}, Id , Va) of the 
diode is equal to — Yc , where Yc is the admittance of the ca\'ity. For a 
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given.setting of the plunger (Fig. 2), F, = Fc(a)) depends only on the fre- 
([Ueney. Thus the frequency and amplitude of oseillation are determined 
in ternis of the direct current by Y{w, Id , l'„) + l%(aj) = 0. When 7j 
is varied both !'„ and w will vary. However, it is not certain whether the 
variation of frequency with Id could be nmde large enough for a practical 
frequency modulation device. Near the optimum frequency and bias 
the frequency will remain relatively constant as Id and Va vary. Thus 
the frequency is primarily determined mechanically by adjusting the 
height of the cavity (Fig. 2) and the amplitude of oscillation is deter- 
mined electrically by the bias. 

4.10 Zener Current 

In sufficiently narrow junctions, where breakdown occurs, the rate 
of generation is an extremely sensitive function of field. So, as in second- 
ary emission, the generation can be highly localized. The diode could be 
made to operate by Zener current rather than multiplication if the p 
region is sufficiently narrow. However, the conditions for negative re- 
sistance would be less favorable. The current /o generated by field emis- 
sion is a function of the peak field Eo . Hence /o and £"0 are in phase, so 
all of the delay has to come from the transit time. For ideal phase rela- 
tions at large amplitudes the bias should be such that the ciUTent is 
generated mainly in a short burst near the voltage peak. Then, if the 
transit time is f of a cycle, the current, le , will flow during the last three 
quarters of the (ycle and power will be dehvered to the ac signal. How- 
ever, the Q and efficiency are considerably lower than for secondary emis- 
sion. The small-signal Q can be varied by varying the bias. A small 
change in bias will change dlo/dEo drastically (since we are on the knee 
of the current-voltage curve). For the frequency 0.75/t, the minimum 
— Q is 25 at small signals and over 100 at large signals. The latter might 
be improved somewhat by increasing the frequency in relation to l/r. 
The efficiency could probably not be raised above about 5 per cent 
without ruining the Q. 

An analysis of the diode operation on Zener current is discussed in 
Appendix E. 
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APPENDIX A 

Derivation of !« 

We first give a simple finalytical derivation and then a physical argu- 
ment. A mort! general proof is given by Shockley (1988). 
Subtract (4) from (3) and use Poisson's equation 

^ = ll« (p - „ + JV) 

ax K 

for ]) — n. Tlic result is* 

^[^£^ + 71=0. 

In other words the (luantity hi brackets is a function of time ordy and is 
the same, at a tjiven time, throughout the length of the diode ^ not 
only in the space-charge region, but also in the ends. The current in the 
ends, or leads, is /, + /«. Let the ends be of sufficiently low resistance 
so that the field can be assumed to vanish there. Then 

f^ + / = /, + /.. 

4t dt 

Averaging ov((r the length W of the space-charge region and using 
C = k/^tW gives 

The physical aigument can be illustrated by Fig. 8(h). Lc^t a sinidl 
pulse of charge 6Q be generated at ;r = 0. It causes a drop in 
field SE = (-iir/KJSQ. If the field at the edges remains constant as the 
pulse moves, the \'olt:ige will drop at a rate vSE = i4Tr/K)vSQ. There- 
fore, if the voltage is (o remain constant, a current 5/„ nuist fiow in the 



* G. Woini'{'icli (lu'lviilc coinmuiiication) has pointed out thitt tlii.'i is a. special 
t-asc of tlie K*'"<'''!il llircc-diineiisional result 



which follows from 



luiil (iiv <^iirl = 0. 



[K aE -. ~| 
-- — -h7 =0 
i-n- dt J 



— K ilE Air -. 
curl //= - — H I 

c dt <; 
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external circuit and increase the voltage by adding charges to the edges 
of the space-charge region. The rate of voltage increase will be hlJC. 
Setting this equal to vhE gives 5/b = SQ/t- = vhQ,IW . 

APPENDIX B 

Instantaneous Pulse 

The equations can be solved exactly for any amplitude and bias in 
the hmiting case of ri ^ 0. As we have seen, the current /o approaches 
an instantaneous pulse. This solution will give reasonable approxima- 
tions to the Q, efficiency and power for actual cases if the effective dura- 
tion of the current pulse is a small fraction of a cycle. We let the fre- 
quency oj be arbitrary but less than 2ir/ci}, so that only one current pulse 
is flowing at a time. The solution is completely specified by the time h 
at which the pulse and the discontinuity m field occur and the magni- 
tude Es of the discontinuity. From {AE} = ~Id/2 and the form of AE{t), 

E, =^{- AE) =^h. (56) 

CO CO 

The condition that h be periodic is 

<^)+|(^^) = (57) 

where E = Eb -\- Va sin cot — Ei{t — ti) in the internal from ts. to 4 -f- 1 
and E = Ef, -{- Va sin coPat all other times. Thus (57) is a quadratic 
equation for Ei as a function of Eb , Va and ti . Only the smaller of the 
two roots is meaningful. When Eb and Va are small compared to unity 
(57) reduces to (37) and is independent of ti . 

To obtain another relation between ii and Es we solve (23) and (25) 
and find the relation between /o and E during the pulse. In the limit of 
an instantaneous pulse, t in equation (23) does not vary during the pulse, 
so dE = — h dt. We can use this to eliminate dt in (25) . Then integrat- 
ipg we have the relation 



hf^] 



/,„..- /o= ^^^!^±^E-^\1+^E\ (58) 



for 7o(0 as a function of E(t) during the pulse. The values of E at the 
beginning and end of the pulse are found by setting /o = in (58). 
Let El be the field at the beginning of the pulse. Then 

El = Eb-\- Va sin coil , (59) 
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and the maximum I,, is 

A,.. = ^^^^'(i + f^,). (00) 

Neglecting teinxs in Ei' in comparison with unity, (58) gives 

;?. = 2fi(i+'^). (m) 

Comliining this with (57) and (CO) we h;ive 

K, + v. sin coil ^ ? /, (l - ^ /,y (02) 

In ])i'actical eases, the term (w?T/6w)/d ?i; Id can be neglected in compari- 
son with unity. Then (03) becomes (55). This, together with the quad- 
ratic (57), determines E6 and h for any Eb and F„ . 
From (00) and (01), 

J {m + 1)^;- 7r"(m + I) j2 /po\ 

Thus the oftective diu'atiou of the pulse M = IV^n^x is 

7r-(m + l)id 

So the pulse becomes sharper as ti decreases and Id mcreases. We can 
estimate the accuracy of the instantaneous pulse approximation by com- 
paring M with the period 27r/a3. 

APPENDIX c 

Iteration Method 

The procedure in detail is the following: Each iteration goes from an 
E{t) through an hit) to a new E{t). In the fust iteration we have to guess 
at E{t). The procedure has been to begin with the E{i) corresponding 
to the instantaneous current pulse and illustrated in Figs. and 9. The 
magnitude ^a of the discontinuity is chosen to satisfy 

(E + fm/2)^'> = 0. 

The time /i when the discontinuity occurs is found from (40). Putting 
this E{t) into (25) and integrating gives a periodic Io(()/hiQ) from which, 
using (23), we find a function F{t) = - AE{f)/h{0). Putting E ^Eb+ Ya 
sin irt - 7(i(0)F(0 into {E + (;m/2)E-) - gives a quadratic equation 
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for /„(()). Only the smallei- of the two voots is menniiigful. From tlie 
known F(t) and Ii.H)} we have a. new E{(). Since it satisfies 
{E + (m/2)E-) - 0, the ln(t)/h{0) calculated horn it in the following 
iteration will be periodic. 

To plot hit) in Fig. 11 we have determined both h(t)/Io{0) and 
/(,(0) from the same E(t). The quantity /o(/)//i,{0) is found from (25) and 
/„(0) is chosen to satisfy (/„) - -2{AE). This relation is automatically 
satisfied if A£'{0 has been determined from /o(/) using (23). However, 
only in the exact solution is it satisfied if Ii,{i) was determined from 
AE(t). To determine /o(0) from AE(I) we have 

\h{0)/ 

where /„(/)//o(0) is found from AE{t) using (25) with E = Et -^ V^ 
sin irt + AE. The value of /a(0) found in this way, in each iteration, 
was compared with the value that makes the following /o(/) periodic. 
When the two values agreed within a specified amount, usually taken 
to be about one per cent, the iteration procedure was terminated. The 
procedure was programmed on an I.B.M. 050 Magnetic Drum Calcu- 
lator. The machine would give. a solution of the required accuracy in 
about four iterations, or forty minutes, on the average. 

In one case, where the current peak occurred at ii ^ 0.92, the itera- 
tion procedure was repeated starting with h -= 0.85. The peak in suc- 
cessive iterations moved from 0.85 to 0.92. In other words, the final 
result was independent of the initial h . However, more itei'atious were 
required for a poor initial choice. 

APPENDIX D 

Effect of I, 

It is convenient to consider the solution for a single cycle extending 
between two current minima. The effect of I, will be important only 
near the ends of the cycle where /„ is small and the term h/h cannot 
1)0 neglected in (52). Let the solution hej(l)ly,(t) whore h{t) is a solution 
for /, = and is the correct solution during most of the cycle. Then 
fit) will be unity during most of the 'cycle and will rise sharply at the 
Ijeginning and end of the cycle. Since the current is always positive, 
each rise in fit) must he less than unity. The condition that the current 
be periodic is A In / + A In h = 0, where A denotes the change during 
a cycle. The change in In h is found by the same procedure that led to 
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(37). If A/ is the total change in / in a cycle we have, instead of (37), 

1„ (I+AA = M!!L±1) r& _ ^.f + 'II vA . (65) 

\2 - A/7 wri L 2 4 J 

Replacing f,)(t) by the correct Kolution /(^)/i)(0 in (52) gives 

at Ti h 

To determine A/ we need to solve this only in the short interval near 
the current minima where / is changing. In practical cases the current 
will have stojiiK-d flowing slightly liefore the end of the negati\'c half of 
the cycle. Hence when /is changing, E + (m/'2)E' will vary a.s Vaut and 






In I '^ I ^ (m+OF.^r ^^.^j 



where 7,„in is the minimum of /o(/) and / ^ is taken at the minimum. 

Substituting (()7) into (Hd) and integrating over the short. intGr\-al where 
/ is changing gives 



The ratio of /„,iu to /,„„x is found by uitegrating (25) from llie (airrent 
minima at K = to the maxinunn at t = ti . At the amplitudes where 
the effect of /,. is important we can use the sharp-iiulsc approximation, 
so E ^ E,. -\- r„ sin w/ during this interval. The quantity /„u,x is given 
by ((33). Thus A/ is found in terms of Ei , !'„ and fi . I'^iuation (1)5) now 
replaces (37) and together with (55) determines /i and Id for a given 
voltage bias and amplitvide. 

APPENDIX E 

Zener Current 

The small-signal case is easily solved. The ac variations hi Ai and En 
are proportional. At large signals it is probably a good approximation 
to say that (a) no ciU'rent is generated until /i'o has risen to a critical 
value, which w(^ shall define as Ec , and (b) thereafter /n will lie such as 
to keep .^11 from rising above Ec . Thus (25) is replaced l)y /n — for 
E <0 and 7? = for h > 0. Equations (23) and (26) remain unchanged. 
The curi'ent /,(/) will flow during the following three intervals: 

(A) Bcgiiming at the time /u when Ea ^ E^ , current liegins to be 
generated and continues to be generated until a lime /j when no more 
carrier si)ace-cliarge is required to keep Eu from rising above Ec , that 
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is, to keep E from becoming positive. The last holes to be generated 
will have been generated before the first have crossed the space-charge 
region. During this interval, h increases from hitn) = to its maximum 
/<;(^i), while /(] jumps at once to its maximum at to and thereafter de- 
clines, reaching zero at ^i . The governing equation is 

r'^- I, = C^ = CV,<^ cos to^. (09) 

dt at 

This can be obtained by physical reasoning from Fig. 8 or by differenti- 
ating (23) and (26). Equation (69) is easily solved, subject to h{io) - 0. 
The time ti and the maximum current Ie(,ii) are found from dle/dt = 
Hit = ti. 

(B) From the time, h , when /« reaches its maximum and h drops to 
zero, until the time to + t, when the first holes to be generated reach 
the i-p^ junction, h remains constant and equal to /e(/i). A constant 
number of holes are moving with constant velocity across the space- 
charge region. 

(C) Between to -{- r and h -\- t the holes generated in the first interval 
are flowing out at the i-p"^ junction. Since the first holes to be generated 
are the first to flow out, /«(/) = Ie(ii) - I,{t - t) during this interval. 

Thus, by solving (69) for hit) in the interval If) ^ t ^ ti we know 
h{t) throughout the cycle. Except in the first interval, E ^ and h = 0. 

The equations have been solved for the frequency 0.75/t or to = 3t/2t. 
The Q is a function of toto , which is determined by the voltage bias and 
amplitude. The Q is negative for 13° < w^o < 90°. The maximum -1/Q 
is about 0.0095 and occurs at about to^o = 40°. The quantity [e/(l -|- e)] 
{Vd/Va) rises almost linearly with wio from zero at 13° to 0.21 at 90°. 
The limitations on Vd and Va would be about the same as for a diode 
operating by multiplication. 
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